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Abstract 

The fluctuation-dissipation theorem is a central theorem in nonequilibrium statistical mechanics 
by which the evolution of velocity fluctuations of the Brownian particle under a fluctuating en- 
vironment is intimately related to its dissipative behavior. This can be illuminated in particular 
by an example of Brownian motion in an ohmic environment where the dissipative effect can be 
accounted for by the flrst-order time derivative of the position. Here we explore the dynamics of the 
Brownian particle coupled to a supraohmic environment by considering the motion of a charged 
particle interacting with the electromagnetic fluctuations at flnite temperature. We also derive 
particle's equation of motion, the Langevin equation, by minimizing the corresponding stochastic 
effective action, which is obtained with the method of Feynman- Vernon influence functional. The 
fluctuation-dissipation theorem is established from flrst principles. The backreaction on the charge 
is known in terms of electromagnetic self- force given by a third-order time derivative of the position, 
leading to the supraohmic dynamics. This self- force can be argued to be insigniflcant throughout 
the evolution when the charge barely moves. The stochastic force arising from the supraohmic 
environment is found to have both positive and negative correlations, and it drives the charge into 
a fluctuating motion. Although positive force correlations give rise to the growth of the velocity 
dispersion initially, its growth slows down when correlation turns negative, and flnally halts, thus 
leading to the saturation of the velocity dispersion. The saturation mechanism in a suparohmic 
environment is found to be distinctly different from that in an ohmic environment. The comparison 
is discussed. 
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I. INTRODUCTION 



One of the fundamental problems in statistical mechanics concerns the microscopic origin 
of dissipation and relaxation of a nonequilibrium system in its course toward equilibrium 
with the thermal bath. It is known that the particle moving in a randomly fluctuating 
medium undergoes Brownian motion. A phenomenological but rather successful approach 
to study the evolution of Brownin motion is by means of the Langevin equation. This 
Langevin equation is a classical equation of motion by phenomenologically adding the terms 
that account for the effects of dissipation and fluctuations under the influence of a fluctuating 
environment. These two effects are ultimately responsible for evolving into thermodynamic 
equilibrium of the Brownian particle with the bath, and are thus related by the fluctuation- 
dissipation theorem (FDT). The FDT is one of the cornerstones of statistical mechanics [1]. 
The static version of the FDT |2j relates the response of a system in equilibrium to action 
of a small external perturbation. The result of it gives a relation between noise correlations 
and susceptibility in frequency domain and the proportionality constant depends on tem- 
perature. The extension of the static FDT to the situation away from equilibrium can be 
realized in terms of the Langevin equation. The Einstein relation, one of the manifestations 
of the dynamic FDT, links the effects between friction and noise correlations in the Langevin 
equation, and thus plays an essential role in stabilizing the dynamics of the Brownian par- 
ticle. The most striking feature of this relation is that it can determine the irreversible 
evolution of the particle from the fluctuations correlation of the heat bath. Thus, the origin 
of this irreversibility of the Brownian dynamics can be attributed to the fluctuating nature 
of the environment. 

A very clear microscopic description, leading to the Langevin equation, within the con- 
text of one-particle quantum mechanics has been presented by Caldeira and Leggett. They 
considered a specific system-environment model that the particle interacts with an envi- 
ronment composed of a infinite number of harmonic oscillators by linear coupling of the 
oscillator and particle coordinates [3]- If the quantum states of harmonic oscillators are 
thermally distributed, the relative importance between quantum and thermal fluctuations 
depends on the temperature under consideration. The effects of environmental degrees of 
freedom on the particle can be investigated with the method of Feynman- Vernon influence 
functional by integrating out environment variables within the context of the closed-time- 
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path formalism [21 HJ O |6l [7] . The more comphcated interaction by considering nonhnear 
couphngs of the particle coordinate can also be studied perturbatively in P|. Under the 
classical approximation where the intrinsic quantum fluctuations of the particle is ignored, 
the Langevin equation can be obtained by minimizing the corresponding effective stochastic 
action. 

At low temperatures, the particle is affected mainly by quantum fluctuations of the 
environment, and it will lead to a new phenomenon, the so-called quantum Brownian motion. 
Many experimentally accessible systems (such as dissipation in quantum tunneling) [3], and 
the problems of quantum measurement theory (such as quantum decoherence of the system 
due to the interaction with its environment) [lOj share the similar dynamics as quantum 
Brownian motion. At high temperatures in which the thermal fluctuations dominate, the 
problem of Brownian motion is described by the classical dynamics introduced in the previous 
paragraph. The relativistic Brownian motion has also been discussed in [9]. 

Here we would like to stress that under the approach of Feynman- Vernon influence func- 
tional, the corresponding FDT can be derived from first principles for a given microscopic 
model in terms of the Green's functions of environment variables. The FDT is a central 
theorem in nonequilibrium statistical mechanics by which the evolution of velocity fluctua- 
tions of the particle under a fluctuating environment is intimately related to its dissipative 
behavior. Thus, in this paper, we wish to explore the dynamics of the Brownian particle in 
the supraohmic environment where backreaction dissipation is governed by the term with 
high-order time derivatives of the position than the one (a flrst-order time derivative) in a 
ohmic case. The known example is the motion of the charged particle under the electromag- 
netic fluctuations at flnite temperature. The non-uniform motion of the charge will emit 
radiation that backreacts on itself through the electromagnetic self-force given by a third- 
order time derivative of the position. The stochastic noise, which encodes the influence of 
quantum/ thermal statistics of the flelds, drives the charge into a fluctuating motion [TTl [T3] . 
It is also known that, in nonrelativistic motion, the dissipation term is signiflcant only on an 
extremely short time scale, which is the time for light to travel across the classical radius of 
the charged particle, ~ lO^^^s for an electron. It is then of interest to flnd out the mecha- 
nism by which the velocity fluctuations saturate, when backreaction dissipation is negligible 
during the evolution, that is, when the particle barely moves in a supraohmic environment. 

The dynamics of a charged particle coupled with the electromagnetic flelds has been stud- 
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ied quantum-mechanically in the system-plus-environment approach. We treat the particle 
as the system of interest, and the degrees of freedom of fields as the environment. The ef- 
fects of fields on the particle is then obtained by integrating out field variables [31 IH El El E] • 
In this approach, the decoherence phenomena of the charged particle under the influence of 
electromagnetic vacuum fluctuations in the presence of the conducting plate has been studied 
in [HI [12]. The evolution of charge's velocity dispersion affected by quantum electromag- 
netic fields near a conducting plate [15] and from the electromagnetic squeezed vacuum [16] 
are also studied. In the latter works, we investigate the possibility of reducing the veloc- 
ity dispersion of the charged particle by tuning the parameters in these quantum states. 
In [IT], the authors have shown that if the particle dose not move significantly such that 
the electromagnetic self-force can be ignored, then the velocity dispersion still reaches a 
constant value at asymptotical times. Thus, constrained by the uncertainty principle, it is 
quite reasonable that the particle cannot extract energy indefinitely from the vacuum state 
of environment fields. It results from the fact that the integral of the correlation function of 
stochastic forces, which will be defined later, over the whole time domain vanishes. Similar 
results are also obtained in Ref. [TSJ. As will be seen later, it will also lead to the saturation 
of velocity fluctuations of a Brwonian particle under a supraohmic environment at finite 
temperature. This work is a follow-up from our earlier investigations [T5l [16 ]. where more 
details on derivation can be found. 

The Lorentz-Heaviside units with h = c = 1 will be adopted unless otherwise noted. The 
metric is r/^'^ = diag(+l, —1, —1, —1). 



II. LANGEVIN EQUATIONS 

The dynamics of a nonrelativistic point particle of charge e interacting with the quantized 
electromagnetic fields can be described by the Lagrangian, 

~1 



L[q,AT] = ^mq2 -^J d^x£y Q{x;q)G{x,y)g{y;ci) + j d^x 



-(9^AT)'+j- At 



in terms of the transverse components of the gauge potential At, and the position q of the 
point charge in the Coulomb gauge. The instantaneous Coulomb Green's function G(x, y) 
satisfies the Gauss's law. The charge and the current densities take the form 

g{x; q(t)) = e6^'\^ - q(t)) , q(t)) = eq(t) ^(^'^(x - q(t)) . (1) 



The density matrix of the particle-field system p{t) evolves unitarily according to 

pitf) = Uitf,U)piU)U-\tf,U) (2) 

with U{tf,ti) the time evolution operator. We assume that the state of the particle-field 
at an initial time is facterizable as p(tj) = Pe(tj) ® pA-riti), and that the electromagnetic 
fields are initially in thermal equilibrium at temperature T = so that its density operator 
takes the form 

PA,(t,,) = e-^^"-/Tr{e-^^-T} , (3) 

where i^Ax is the Hamiltonian of the free fields. 

After integrating out the degrees of freedom of the fields, the Langevin equation is ob- 
tained [IH] . 



where 

P(f\ = -hp ( S^' 

dt 

with the noise-noise correlation functions, 



m = -ne f^^'^-^'WvM e'W (5) 



= 0, (r(t)e(O) = \G%ut),c^{t')-t-t'] , (6) 



and 



hG%{x-x') = ie(t-t'){[A'T.{x),A{{x')]) , (7) 
hG%{x-x') = i({A^,(x),AUx')}>, (8) 

are the Green's function of the electromagnetic potentials at finite temperature. It is seen 
that the infiuence of the electromagnetic fields are expressed by an integral of the dissipation 
kernel over the past history of charge's motion, and by a stochastic noise ^ that drives the 
charge into a fiuctuating motion. As it stands, Eq. Q is a nonlinear Langevin equation with 
non-Markovian backreaction, and the noise depends in a complicated way on the charge's 
trajectory because the noise correlation function itself is a functional of the trajectory. 

The fluctuation and dissipation effects of the electromagnetic fields on the motion of the 
charged particle are associated with the kernels and G^ respectively. They in turn are 
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linked by the fluctuation-dissipation relation where the Fourier transform of the fluctuation 
kernel is related to the imaginary part of the retarded kernel as follows 



G'^[q(t),q(t');cu] = Im{G^i[q(t),q(t');u;]}coth 



(9) 



This relation is established from first principles. The explicit expression of the kernels 
and G^ will be introduced when they are needed later. 



III. VELOCITY FLUCTUATIONS 



The nonlinear, non-Markovian Langevin equations are far too complicated to proceed 
further without any approximation. The appropriate approximation for nonrelativistic mo- 
tion is the dipole approximation, which amounts to considering the backreaction solely from 
the electric fields. The Langevin equation under the dipole approximation reduces to 



from Eq. (|4]) and 



mq\t) + e' / dt' gl{t - t') q\t') = fl{t) 
Jo 



fl{t) = -heC{t). 



(10) 



(11) 



The retarded Green's function in the dipole approximation is denoted by qr, and can be 
expressed in terms of the spectral density p as 

""^ dk 



p'^{k) sin(A;r) . 



(12) 



In the isotropic thermal bath, the spectral density takes a simple form 

k 



(13) 



The accompanying noise-noise correlation functions due to the fluctuations of the electric 
fields at finite temperature are given by 

92 



with 



Ul{t)n{t'))=he^—g^^,{t-t') 



dk 
2^ 



k'^p'^{k) coth 



dtdt 
/3hk 



cos(A;i 



(14) 



(15) 
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Again the noise kernel g^^ can be seen related to the dissipation kernel g^^ via a fluctuation- 
dissipation relation under the dipole approximation, derived from Eq. ([9]). After carrying 



out the integration in Eq. (10) the Langevin equation becomes physically more transparent, 

mA\t)-^V{t) = rs{t). (16) 

The non-uniform motion of the charge results in radiation that backreacts on the charge 
itself through the electromagnetic self-force. This backreaction occurs at the moment when 
radiation is emitted [131 El US]- Thus, it may lead to short-distance divergence in the 
coincidence limit due to the assumption of a point-like particle. This ultraviolet divergence 
must be regularized to have a finite and unambiguous result. The divergent part is absorbed 
by particle mass renormalization, = m + e^A/Svr^, where A is the energy cutoff scale 
related to the inverse of the charge's wavepacket width. It essentially quantifies the intrinsic 
uncertainty of the charged particle. Then the finite backreaction effect is given by the 
well-known result, a third-order time derivative of the position |20j . 

For a nonrelativistic particle, the introduced cutoff wavelength should be much larger 
than its Compton wavelength, namely A^^ ^ \c = h/m^c^, and in turn much greater 
than the classical radius of the charged particle, r^. = /m^c. Thus, when the time scale 
associated with the cutoff frequency, t ~ 27r/A, is much longer than the characteristic time 
scale Te ~ Te/c, an extremely short time scale, the electromagnetic self-force can be safely 
ignored. If we assume that the charged particle starts off from the rest at t = 0, then the 



solution to the Langevin equation (16) has a very simple form 



vm = — fdsfM, (17) 



and the corresponding velocity dispersion is given by 



{Avm) = —J ds [ ds' {fi{s)f:{s')) . (18) 

Jo Jo 

We see that the evolution of velocity fluctuations is governed by the force-force correction 
function, whose finite temperature contribution (/1(t)/*(0))/3 is given by 

Cff(t) = (/:(r)/:(0))^ 
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Cff(t) 




FIG. 1: The correlation function Cff{'t) normalized with is drawn as a function of r in 



units of 



where its vacuum part has been subtracted, ({n, z) is the nth derivative of the zeta function. 
The correlation function Cff{t) is plotted against the time difference t = t — t', scaled by 
(3h in Fig. [Tj It shows that the stochastic noise tends to have positive correlation for small r, 
and then Cff{t) turns negative at r ~ 0{f3). It means that if we measure the electric field 
at certain time and find out the field points to one particular direction, then the moment r 
later, within the time difference of order 0{(3), if we perform a similar measurement again 
at the same location, there is high probability that the direction of the electric field remains 
the same. However, when the time difference between measurements is greater than 0{l3), 
we probably have the result that the electric field points to the opposite direction instead of 
the same direction. Another way of understanding its significance is to compute the mean 
power Psit) done by the stochastic force fs(t), 



m 



ds {mfs{s))p 



r Jo 



m 



1 he^ 



Im|^/'(2)(1 + 



ds Cppit — s) 



r Jo 



(20) 
(21) 



where il)'^"'\z) is the rath derivative of the digamma function il){z). Again here only the 
finite temperature contribution is considered. We see that the mean power is positive at 
early times where the stochastic force does positive work to the charge such that its velocity 
dispersion or kinetic energy increases with time. The power reaches the maximum value at 
the time when correlation turns negative. Later on, the negative correlation makes the work 
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p. it) 




(At.?{t)) 



4 ph 



FIG. 2: The time evolution of Ps{t) and {Avf{t))f3 is plotted where the time t is in units of 
These quantities are normalized with 



I3h- 



3«i% respectively 



done by the force less positive and then the mean power Psit) eventually approaches zero, 

^TT^ t t 



Then the finite temperature part of the velocity dispersion {Avf{t))i3 can be obtained by 

2 



m 



ds P Js) 



T Jo 



he" 



D t 



t:^ ph. 



In Fig. [2jwe show the time evolution of Ps(t) and {Avf(t))p respectively. It is seen that 
the averaged power Ps{t) given by stochastic forces on the charged particle is positive for 
all times. This implies that the thermal bath keeps pumping energy to the particle during 
its evolution and increases the velocity dispersion. The presence of the negative correlation 
slows down the rate of energy transfer, and finally leads to the vanishing Pg at asymptotic 
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times where velocity fluctuations are saturated. It implies that the integral of Cff{t) over 
the whole time domain vanishes, i.e. 

poo poo 

/ dr Cppir) = 0= dr Cff{t) . (22) 

Jo J -oo 

Thus, the contribution from the positive correlation for the small r regime exactly cancels 
that from the negative correlation for the larger r regime. It leads to the saturation of the 
velocity dispersion even though the self force is insignificant in this supraohmic case. 



To understand Eq. (22) and its consequence from a different aspect, we may rewrite the 
correlation function in terms of its Fourier transform 

dk / 1 \ 

"'^ '2 Jin \ I ^ \ -iki 



Then the integration of the correlation function over the whole time domain becomes 

,2 ii/- ^ ' 1 



dr Cff{t) = -j^ dk ep'\k) y-^f^j—^j m 

If kp^^ik) behaves like k^~^ with / — 1 > 0, then Eq. (22) holds. Thus, whether the integral of 
the correlation function over the whole time domain vanishes or not relies on the behavior 
of kp^^{k) in the zero momentum limit, which in turn depends on the spacetime dimension 
of the system, the dispersion relation of the environment field, and the coupling between 
the particle and environment. 

In electrodynamics, the dynamics of a charged particle is governed by a local, gauge 
invariant interaction with the electromagnetic fields. Under the prescription of minimal 
coupling, the transverse component of the vector potential couples to the charged current 
density, which is proportional to the time derivative of particle's position. This derivative 
coupling gives rise to the fact that the electromagnetic self-force should be given by a third- 
order time derivative of the position, and it is called supraohmic dynamics [31 [7j. From 



Eq.(13), we see that / = 3 and hence Eq. (22) holds. In this supraohmic case, the saturation 
of velocity dispersion can be achieved due to negative force-force correlations even when the 
backreaction effect of the self-force is found negligible in the evolution. 

This is in striking contrast with the Brownian motion under an ohmic environment, 
characterized by the dissipative backreaction of the first-order time derivative of particle's 
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position. It involves the coordinate coupling of the particle with the environment [3l[7], and 
then / = 1. The positive force- force correlation pQ has been shown to drive the growth of the 
velocity dispersion linearly in time within a time scale shorter than the relaxation time. It is 
also found j21] that the integration on the force-force correlation function for the whole time 
regime then has a nonzero value with Z = 1, so the energy transferred from the environment 
to the particle during the whole evolution of particle's motion never slows down to halt. Then 
the dissipative backreaction must be taken into account to counterbalance the effect from the 
force fluctuations in order to finally stabilize the value of the velocity dispersion. Moreover, 
from the aspect of energy balancing between dissipative and fluctuation backreations, the 
stronger dissipation is expected to occur in the subohmic case with Z < 1, because the energy 
transfer rate increases even faster at late times [7]. Thus, the mechanisms to stabilize the 
velocity dispersion of a particle in a fluctuating subohmic, ohmic, or supraohmic environment 
are rather different. 



IV. SUMMARY AND CONCLUDING REMARKS 

In this paper, the evolution of velocity fluctuations of the Brownian particle in a 
supraohmic environment is studied by considering the motion of the charged particle coupled 
to electromagnetic fluctuations at finite temperature. The Langevin equation of the parti- 
cle incorporates the effects of fluctuation and dissipation backreaction in a self-consistent 
manner. In particular, the backreaction in a form of the electromagnetic self-force on the 
charge is a third-order time derivative of the position in this supraohmic environment. On 
the other hand, the thermal fluctuations of the electromagnetic fields are manifested as 
stochastic noise. Its correlation tends to be positive at shorter time difference r, and then 
turns negative at r ~ 0{(3). We show that the integration of the force-force correlation 
function over the whole time regime vanishes. Throughout the evolution, the self-force can 
be argued to be insignificant for the charge having no significant motion. Then, the positive 
correlation contributes to the growth of the velocity dispersion initially, its growth slows 
down when correlation becomes negative, and finally halts where the velocity dispersion 
reaches a constant at asymptotical times. It is a rather different saturation mechanism for 
the velocity dispersion of the Brownian particle in an ohmic environment. 

Additionally, the saturated value of the kinetic energy of the charge due to the thermal 
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electromagnetic fluctuations at temperature T — (5 ^ can be estimated by 

lm,(A.?(oo)), ~ « ^ ~ 10- (*f ) ' keV , (24) 

with the flne structure constant a = /he. The a dependence imphes that the saturated 
value of the velocity dispersion in this supraohmic case relies on the coupling between the 
system and the environment. The issue on how the velocity dispersion of the Brownian 
particle is stabilized in the subohmic environment deserves further investigation. 

Thus, it comes to no surprise that the above saturated value of the velocity dispersion 
roughly becomes the value given by the equipartition theorem for a thermodynamic system 
as long as /c^T ~ rrirC^ /a ~ ficT~^ because in that limit, the temperature is the highest 
energy scale in the system. However, in order to correctly describe the system in such a limit, 
a relativistic/ficld-theoretic formalism is needed for appropriately describing the dynamics 
of the charged particle. This is beyond the scope of the current investigation. 
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